Noise spectroscopy using correlations of single-shot qubit readout 
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A better understanding of the noise causing qubit decoherence is crucial for improving qubit 
performance. The noise spectrum affecting the qubit may be extracted measuring dephasing under 
the application of pulse sequences but requires accurate qubit control and sufficiently long relaxation 
times, which are not always available. Here we describe an alternative method to extract the 
spectrum from correlations of single-shot measurement outcomes of successive free induction decay 
measurements. This method only requires qubit initialization and readout with a moderate fidelity 
and also allows independent tuning of both the overall sensitivity and the frequency region over 
which it is sensitive. Thus is is possible to maintain a good detection contrast over a very wide 
frequency range. We discuss using our method for measuring both 1/f noise and the fluctuation 
spectrum of the nuclear bath of GaAs spin qubits. 
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Understanding and reducing decoherence of qubits is of 
great interest for quantum computing. Pulse sequences 
like Hahn's spin echo [l[ and more advanced sequences 
such as CPMG 0, concatenated (CDD) or Uhrig dy- 
namical decoupling (UDD) Q have been demonstrated to 
decouple a qubit from its noisy environment and there- 
fore reduce decoherence. In addition to enhancing de- 
phasing times, such pulse sequences can also be used for 
noise spectroscopy by observing the resulting dependence 
of qubit coherence on the evolution time. Ref. [5| shows 
that for Gaussian noise, the spectrum's moments can be 
obtained using UDD. Yuge et al. proposed a sequence 
of equidistant 7r-pulses to reconstruct the dephasing noise 
spectrum. A different approach relics on separating the 
times scales of the noise correlation and extracting long- 
time correlations via direct measurement while investi- 
gating short-time correlations with pulse sequences 0. 
However, pulse-sequence based noise spectroscopy is sub- 
ject to certain limitations: For fixed pulse sequences, the 
frequency region over which this technique is sensitive 
is proportional to the inverse evolution time, but longer 
durations also increase the total decoherence. This rela- 
tion makes it hard to probe low frequency noise: by the 
time the frequency region of interest is accessible, the 
qubit is fully dephased leaving no measurement contrast. 
Depending on the details of the spectrum, this problem 
may be to some extend circumvented by adding more 
pulses and appropriate choice of their timing. However, 
this strategy will eventually be limited by 7r-pulse er- 
rors and the T\ time of the qubit. Consider for example 
experiments probing the flux noise in superconducting 
qubits. Direct measurement of the state allows for de- 
tection of low-frequency spectral content < 1 Hz 0,0 
while CPMG and UDD sequences have been used to re- 
veal high-frequency content > 0.1 MHz [l0 thus leaving a 
gap inaccessible with those methods. Similar limitations 
are encountered in GaAs electron spin-qubits, where the 
performance of decoupling sequences was likely limited 



by imperfect control pulses rather than intrinsic noise 
caused by nuclear spins 11 1. 

In this Letter we propose an alternative method for 
determining the noise spectrum at intermediate frequen- 
cies without requiring 7r-pulses, and which is capable of 
probing the noise spectrum over an extremely wide fre- 
quency range. It is based on correlating (near) single- 
shot free induction decay (FID) measurements. A typical 
measurement cycle is depicted in Fig. [1] and consists of 
initializing a a x eigenstate, performing a free evolution 
under the influence of the noise process for time t, a pro- 
jective measurement of the final state, and repeating this 
process after a delay time At. Averaging over many such 
measurements allows computing the correlation between 
consecutive measurements as a function of At. A similar 
method was employed in the experiments of Rcilly et al. 
[l3 |. where correlations of temporal averages instead of 
single-shot readout were used to determine the nuclear 
spin noise spectrum seen by an electron spin qubit up to 
1 kHz. A different approach uses dynamical decoupling 
techniques to introduce an effective delay time and mea- 
sure correlations of the spin bath of a nitrogen vacancy 
qubit [13 1 . 
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FIG. 1. (color online) Measurement cycle: Qubit initialized 
(I), evolution (E) for time n and measurement (M) of the 
outcome. The delay time between two evolutions is At. While 
the initialization and measurement time set a lower limit on 
At, the idle time can be varied to adjust the delay between 
evolutions times. 
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After discussing the relationship between noise spec- 
trum and single-shot correlations in general terms, we 
consider two specific examples from recent experiments 
to demonstrate the usefulness of our method. 

The energy relaxation time of our qubits is assumed 
to be much longer than the evolution time of interest 
and we therefore only consider pure dephasing. This as- 
sumption is adequate for many experiments on electron 
spins [l2, 14 and for superconducting qubits at op- 
erating points where they are linearly sensitive to flux 



noise 



lol 17 -Til- The qubit Hamiltonian therefore is 



H 



[n + p(t)] 



(i) 



where f2 is the qubit energy splitting and j3{t) represents 
a classical noise process with (/3) = 0, which we assume to 
be Gaussian. We first set O = and will discuss the gen- 
eral case later on. The initialized state is a a x eigenstate. 
A subsequent measurement of a x has an expectation 

value P = cos[A$(t)], where A$(t) = J^/jf PW) dt' 
is the phase accumulated during time r due to noise. 
The autocorrelation function of the single-shot outcomes 
is then given by 

(P / (t).P / (t + At)) = I e -^(A*(t) + A$( t +A t )) 2 > 

+ I e -i((A*(t)-A<I>(t+At)) 2 } ( 2 ) 

= 2 e 2 + 2 e ' 

and takes values between 1 (perfectly correlated) and 
(completely independent). For simplicity, we considered 
symmetric measurement cycles with n = T2 in writing 
Eq. [2] Asymmetric evolution times can be easily intro- 
duced to Eq. [3] and lead to reduction of the autocorrela- 
tion function due to uninteresting FID-typc effects. We 
express the exponents x± m terms of the noise power 



spectrum Sp(u>) = 



j At 



(/3(t)/3(t + At))dAt via 



(A$(f)A$(f + At)) = / — Sg(u>) 
Jo n 



F(lut) 



where 



F(utTi) = 4cos(wAf) 



sin 



V 2 



(3) 



(4) 



is the filter function analogous to those of de- 
coupling sequences @. The variances x± = 
((A$(i) ± A$(f + At)) 2 ) entering Eq. d] are calculated 
straightforwardly. For r >• T£ > the terms ((A$(f)) 2 ) are 
much larger than unity, but \- (x+) can be small be- 
cause of correlations (anticorrelations) leading to a par- 
tial cancellation of the variances by the (A$(t)A$(f + 
At)) term. In this case, only one of %_ and x+ can be of 
order unity, so that the other one is negligible because its 



exponential vanishes. Since correlations are more com- 
mon than anticorrelations, the dominant term will typ- 
ically be X- an d thus can be measured directly. This 
behavior can also be understood from the filter function: 
X+ has a filter function with larger weight at low frequen- 
cies where the spectrum is typically largest, whereas the 
filter function for x~ vanishes quadratically for w — > 0. 

For At = t, the filter function of x~ is identical to that 
of a spin echo experiment where F(ujt) = 8 sin 4 (2cjt/4) 
and the effect of the reinitialization on x~ is thus equiv- 
alent to a 7r-pulse. Note however that to reach this limit, 
initialization and readout need to be negligibly short. 
While we have assumed fl = so far, a nonzero adds 
oscillations oc cos (2f2r) to the first term of Eq. O How- 
ever, the second term containing X- remains unaffected 
for T\ = T2, so that our results remain valid for f2 7^ 0. By 
varying the delay time, At, we can shift the weight of the 
X--filter function in the frequency domain and thus tune 
the highest sensitivity to the frequency range of interest. 
Note that for wt -c 2, the evolution time r controls the 
gain factor and thus the overall sensitivity. The ability to 
adjust both the sensitivity's overall magnitude and posi- 
tion mark the strength of our method: suitable choice of 
evolution and delay time allows for maintaining a good 
measurement contrast over a wide frequency range and 
thus enable the investigation of spectral content. To max- 
imize the range of applicability of our technique, fast ini- 
tialization and measurement would be beneficial as they 
set the lower bound on At. However, low single shot fi- 
delities can be averaged out and slow initialization can be 
replaced by projective (non-demolition) measurements. 

As a concrete application of our method, we consider 
electron spin qubits based on gate-defined GaAs quan- 
tum dots. For two-electron spin qubits in double quan- 
tum dots, the requirements of fast initialization and read- 
out have already been shown [l5|, [20[, and they could 
likely be achieved with similar techniques for single-spin 
qubits. For these qubits, the hyperfine interaction cou- 
ples the electron spins to ~ 10 6 nuclei in their respective 
dots. The resulting Overhauser field B nuc contributes to 
the Zeeman splitting of the electrons and its temporal 
fluctuations are the main source of dephasing. For two- 
electron singlet-triplet qubits, the difference of the Over- 
hauser field in the left and right dot AB = £?n UC — B^ uc 
acts in the same way [HI, [2l| . 

For this qubit system, /3(f) = ^-AB z (t), where g = 
—0.44 is the effective gyromagnetic factor for electrons 
in GaAs and fiB is Bohr's magneton. The temporal fluc- 
tuations in the dots are due to spin diffusion processes 
that enable distant nuclear spins to exchange polariza- 
tion. Treating these processes as classical diffusion, one 
finds that the noise spectrum of the Overhauser fields is 
proportional to lu~ 2 . Theory 22j and experiments in- 
dicate an exp(— t 4 ) decay, which implies that the w -2 
dependence is limited at approximately 10 kHz. How- 
ever, available data (ill |23| is insufficient to quantita- 
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tively characterize the high frequency behavior. Ref. [24[ 
argued that the nuclear spin coupling strength sets an up- 
per cutoff to this behavior and modeled the rolloff as a 
reduction factor of the form exp [— (lu/loe) 7 ]- To demon- 
strate that our method could verify this cutoff behavior, 
we consider a spectral density of the form 



S 



1 



2 x exp 



LJ E 



"< 1 



(5) 



The lower cutoff frequency (u>l) coarsely mimics the low 
frequency behavior of the diffusive dynamics. It ensures 
a finite r.m.s. Overhauser field {y/ (ABf)) typically on 
the order of a few mT (25[. The cutoff frequencies, lot, /2tt 
and loe/^tt, are set to 0.1 Hz and 10 kHz respectively [24[ 
and Sq determines the overall noise level. 

The autocorrelation function of Sab z is shown in Fig. [2] 
(a). The computed correlation functions show the ex- 
pected behavior: for short delay times, the probability 
that the two measurements return to the initial (singlet) 
state are perfectly correlated. By increasing At, this cor- 
relation decays to zero where consecutive measurement 
outcomes are completely independent. For short evolu- 
tion times a larger change in ft is required for decorrela- 
tion, so that correlations persist to longer times. 

It is useful to consider analytical approximations for 
the relevant term of Eq. [2]). We assume that the 
evolution time is short compared with the relevant fre- 
quencies so that tlue *C 1. With this approximation we 
obtain 



1 



-£ C W 2 r 2 S At 2 
1 - 2 *> 2 - 2 SoAt 



At < Ug 1 , 



where 



^r- is a constant and a = V 



J 1 < At < Uj}, 
-c 2 ^-((AB 2 )) 2 At»Lu Z \ 

(6) 

h it, ^uiiD»»ii» dii^ u, — x ^ 7 f lj is given 

in terms of the Gamma function T(x). We find a = 1 
(y/Tr/2) for 7 = 1 (2). Fig. [2] (b) compares these approx- 
imations (dashed, dash-dotted and dotted line) with the 
exact, numerical result (solid line). 

We note that in the regime of short evolution times and 
small noise variations r 2 ((/3(t)/3(t + At)) - (/3 2 )) < 1, 
the single-shot correlation depends linearly on the noise 
autocorrelation 12J, which provides direct access to the 
noise spectrum. We find that 



(P(t)P(t + At)) 



1 



. {m p(t + At))--(p 2 ). (7) 



In order to use our method for noise spectroscopy and, 
more specifically, to investigate the predicted cutoff be- 
havior, we wish to identify the exponential's power law 
7. For constant r, a difference between the autocorre- 
lations for exponential and Gaussian cutoff appears for 
delay times on the order of the inverse cutoff frequency. 
This frequency also corresponds to the region where the 
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FIG. 2. (color online) (a) Autocorrelation of return proba- 
bility P for evolution times from 50 ns (black) to 5 fis (blue) 
and exponential cutoff (7 = 1) for the spectrum in Eq. [5] 
Note that for the dotted part of the blue curve At < r and 
thus the correlation function does not exist, (b) Exponent of 
the dominant term of Eq. [2] Note the three different regimes 
At < 1/u>e, I/Ve < At < 1/ujl and l/u) L < At corre- 
sponding to At 2 , At 1 and At behavior, respectively, as dis- 
cussed in the text. Short evolution times tlue -C 1 result in 
a prefactor r 2 and show a qualitatively similar behavior. 



transition between At 2 and At-behavior takes place. For 
long delay times, the two curves merge as in this region 
Eq. [6] does not depend on the exact form of the spec- 
trum, but only on the r.m.s. value of the Overhauser 
field. The delay time needed to distinguish different cut- 
off behaviors is thus on the scale of 1/uje- Note that in 
this regime no difference between the logarithm of Eq. [5] 
and — \x~ can be detected which retrospectively justifies 
the assumption that the autocorrelation is dominated by 
X-- However, for a constant evolution time, X- varies 
over a large range as a function of At, thus limiting the 
range of delay times for which the autocorrelation is mea- 
surable. We therefore compensate for changes in x- due 
to a varying At by simultaneously adjusting the evolution 
time r. In Fig. [3] (a) we chose r so that the linear term of 
Eq. [6] always equals 2. The resulting curves allow a clear 
distinction between different cutoff exponents 7, while a 
change of the cutoff frequency w e has a qualitatively dif- 
ferent effect. Note that the signal level of ~ 0.1 allows 
a good measurement contrast compared to realistic noise 
levels over a large range of At. 

We now consider using our method for characterization 
of 1//" noise which has been found to be an important 



4 




0.0 



^^^^ no cutoff 
7=1, 00 /271=10 kHz 

E 

7=2, CO c /2Jt=10 kHz 

1 E 



7=1, co /27i= 5 kHz 

E 

7=2, m /2tc= 5 kHz 

1 E 



250 



500 10 
At (us) 



10 




FIG. 3. (color online) (a) Autocorrelation function for the 
Singlet-Triplet qubit. Adjusting the evolution time allows for 
detection of the crossover and slope of curves with different 
cutoff behavior, thus providing well distinguishable charac- 
teristics on an experimentally measurable scale. The blue 
curve shows the autocorrelation in the absence of a cutoff. 
For At > 600 /is the scale is semilogarithmic. Note that the 
plateau at large At directly reflects the prefactor of the 1/uj 2 
region of the spectrum, (b) Autocorrelation function for the 
l/f a spectrum. For r <C At <C exact 1// noise leads 

to a constant result. Deviations from this behavior (red and 
blue curves) allow for identification of a. 



source of dephasing for several types of superconducting 
qubits @, [IE E3, 111- For a fixed r < At and a > 1, 
we obtain a decay of the form \- ^ A/( Q_1 ) while for 
a = 1, X- ln(A//r)r 2 . As a is reduced below 1, the 
At dependence of x~ becomes increasingly weak, so that 
a crossover to a constant x~ corresponding to white noise 
is obtained. By adjusting r oc Aiexp [W (— ^fr)] , where 
W{x) is the solution to x = wcxp(w) and c is a constant 
setting the overall noise level, exact l/J-noise leads to 
a constant result for r <C At, so that small deviations 
from 1// in Fig. |3] (b) can be detected with high sen- 
sitivity over a very large frequency range. This allows 
in principle to probe arbitrarily low uj compared to the 
lower bound of ~ 0.1 MHz of Ref. which was lim- 
ited by the T2 time of the qubit. For superconducting 
qubit systems, initialization is realized trough relaxation 
making our previous assumption of fast initialization in- 
valid as the initialization time is on the order of T\. To 



overcome this problem, initialization can be replaced by 
measurement which is fast and also provides the infor- 
mation needed to evaluate the autocorrelator. 

In conclusion, we have developed a method for noise 
spectroscopy that fills the frequency gap remaining 
between direct and pulse sequence-based noise spec- 
troscopy. Moreover, it enables the identification of dif- 
ferent cutoff behavior and frequency-dependencies in spin 
and superconducting qubits, respectively. As FID pulses 
and nanosecond-control of qubits are well established 
15 . 26-2^|, experimental implementation should be fea- 
sible. Introducing additional operations between the two 
evolution times might reveal whether the qubit state has 
non-negligible effects on the noise bath resulting in dif- 
ferent autocorrelation functions. 
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